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1. Introduction 

The Square Kilometre Array is a faeility whieh will be able to provide huge advanees in several 
areas of astronomy and eosmology. It will be able to map the large seale distribution of galaxies as 
well as measure the weak gravitational lensing signal from far away objeets. The information from 
these probes is not fully independent as both the LSS and WGL signals depend on the large seale 
distribution of matter in our Universe. 

The eombination of the information eneoded in the large seale distribution of galaxies and in 
the signal from weak gravitational lensing ean signifieantly inerease our knowledge of the Universe 
we live in. Weak Gravitational Lensing (WGL) does not respond to the bias of galaxies and henee 
the eombination of the signal from WGL and LSS ean help eonstrain the bias of galaxies traeing 
the Large Seale Strueture. Furthermore WGL usually requires photometrie redshifts to yield an 
estimate of the redshift distribution of sueh galaxies. The eross eorrelation between WGL samples 
with a speetroseopie redshifts samples ean help ealibrate the redshift distribution of WGL galaxies. 
Finally if we eonsider theories of modified gravity, WGL eombined with the Redshift Spaee Dis¬ 
tortion (RSD) signal in the LSS survey ean together probe potential modifieations to the Poisson 
equation as well as the equation whieh determines how light bends in the presenee of matter in 
separate ways. This breaks a very powerful degeneraey in the common parameterisations of devia¬ 
tions from GR, hence their joint constraints are orders of magnitude more powerful than individual 
constraints. 

In this chapter we take surveys equivalent to those possible with SKA phase 1 as well as SKA 
phase 2 and estimate how well the combination of probes can calibrate photometric redshifts and 
improve the figures of merif for Dark Energy and Modified Gravify. We consider several scenarios, 
nofably we describe as SKA 1 early scenarios which are surveys which will be available during 
a consfrucfion phase of fhe SKA before phase 1 is complefed. We assume fhaf fhere is bofh a 
confinuum and a line survey for fhe Weak Lensing and fhe Large Scale Sfrucfure respecfively. 
The Confinuum survey would have longer baselines and would be focused on weak lensing and 
is similar fo fhe survey ouflines in fhe Weak Lensing chapfer in fhis science book. The redshifls 
for fhese galaxies would come from mafching fo optical galaxies wifh phofomefric redshifls. The 
LSS surveys are assumed fo be galaxies found in line emission mainly from fhe signal presenf in 
fhe core of fhe SKA. We make lilfle disfincfion befween fhe fechnologies needed for such surveys 
buf note fhaf a proper UV disfribufion is needed for a Weak lensing survey and some necessary 
sensifivify is needed for fhe line survey. These SKAl early surveys are effeclively 1000 and 5000 
sq degree surveys assuming fwo fhirds of fhe SKA phase 1 sensifivify. We also assume fwo SKA 
phase 1 surveys which would reach signals of around 100 /tJy in fhe case of a line survey. The 
Weak lensing survey is assumed nol fo be larger fhan 5000 sq degrees in fhe case of phase 1 as 
anyfhing more would be infeasible given fhe sensifivify of fhe inslrumenl. We assume a larger LSS 
survey of 30000 sq degs wifh phase 1 obfaining galaxies and redshifls. For phase 2 we assume bofh 
WGL and LSS surveys covering fhe same amounf of fhe sky wifh 5000 and 30000 square degrees 
respectively where increased area Irades off wifh decreased deplh. 

These surveys are summarised in Table 1 and 2 and Iheir respective redshifl disfribufion is fhe 
same assumed in olher chapters of fhis science book. 
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2. Formalism for Cross-Correlations 


In the following we will provide forecasts for various cross-correlations of cosmological 
probes using the Fisher Matrix formalism (see Heavens (2009) for an overview). The Fisher matrix 
is defined as 



( 2 . 1 ) 


where L is the likelihood and where pa are cosmological parameters. Assuming a Gaussian like¬ 
lihood for the data and a cosmology-independent data covariance, the Fisher matrix is given by 
Tegmark, Taylor & Heavens (1997) 
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For any unbiased estimator the Fisher matrix provides a lower bound on the marginalised error 
of a parameter pa, via the Cramer-Rao inequality, Apa > s/{F^^)aa- The data vector, 
consists of angular power spectra, (^) ^ function of multipole i, for a given combination of 

probes, X and Y, and a pair of redshift bins i and j. The covariance matrix in Eq. (2.2) is assumed 
to be Gaussian and takes into account shot noise as well as cosmic variance contributions, see e.g. 
Takada & Jain (2004) for the weak leasing case. 

Following the notations of Joachimi & Bridle (2010), the angular power spectra of the cross¬ 
correlations between the various quantities related to the gravitational potential or the matter den¬ 
sity can be written in a generic way using the Limber approximation. 
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where x is comoving distance, and Xhor the comoving horizon. The matter power spectrum is 
denoted by Pg and the comoving angular diameter distance by fK{x)- The kernels for the different 
cosmological probes are given by the following equations. 
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(2.5) 

(2.6) 


Here, a denotes the scale factor, ph) the redshift probability distribution for a galaxy sample i 
(either a broad tomographic bin or a narrow range defined via spectroscopic redshift), and bg the 
galaxy bias which can vary as a function of scale and redshift. CMB leasing has a single source 
distance, x*^ to the last scattering surface. 

Note that the Limber approximation assumes that the kernels involved are broad in the line-of- 
sight direction, which breaks down in the case of spectroscopic clustering information. In this case 
we replace the formalism outlined above with the exact calculation as detailed in Padmanabhan et 
al. (2007). 
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Survey 

Area [deg^] 
Photo-z 

Area [deg^] 
Spec-z 

rig [arcmin 
Photo-z 

Ug 

[arcmin 

Spec-z 

DE EoM inc. Planck 
WGE ESS WGExESS 

SKAl early 

1,000 

1,000 

4.8 



- 

3.0 

- 

- 

SKAl early 

5,000 

5,000 

1.6 



- 

3.6 

- 

- 

SKAl 

5,000 

5,000 

4.2 



1.8 

16 

0.05 

19.2 

SKAl 

5,000 

30,000 

4.2 



4.7 

16 

1.8 

956 

SKA2 

5,000 

5,000 

44 



22 

182 

320 

6362 

SKA2 

30,000 

30,000 

20 



10 

425 

1042 

13150 


Survey 

Area [deg"^] 
Photo-z 

Area [deg"^] 
Spec-z 

MG EoM inc. Planck 
WGE ESS WGExLSS 




SKAl early 

1,000 

1,000 


1.9 

- 

- 




SKAl early 

5,000 

5,000 

2.0 

- 

- 




SKAl 

5,000 

5,000 

8.0 

4.7 415 




SKA2 

5,000 

5,000 


69 

257 7429 




SKA2 

30,000 

30,000 

118 

583 29947 




Table 1: (Top) Summary of SKA survey forecasts and Dark Energy Figure of Merit results including Planck 
priors.(Bottom) Summary of survey forecasts For the modified gravity Figure of Merit from SKA phases 
one and two. The MG FOM here is defined in the same way as the DE FOM with the parameters of interest 
being the constants modifying the relation between the potentials. 


We forecast constraints for our weak leasing (WGL) surveys alone, {Cee^(£)}, large-scale 
structure (LSS) surveys alone (galaxy clustering including redshift-space distortions), {CnJ\l)}, 
and the joint WGLxLSS analysis including cross-correlations, {CgP {£), Cne (^) j Cm^ (^)} ■ We also 
provide an assessment of the constraining power of CMB lensing with SKA LSS, 
and CMB lensing with SKA WGL probes, 


3. Weak lensing - galaxy position cross correlations 

The SKA survey will provide very competitive measurements of multiple cosmological probes. 
Combination of different probes allows the breaking of degeneracies between parameters which 
provides better control of systematics and better constraints on parameters of cosmological interest 
than any individual probe can achieve. 

In this section we consider the combination of an SKA Weak Lensing (WGL) survey (with 
photometric redshifts provided externally) and an SKA galaxy position i.e. Large-Scale Structure 
(LSS) survey with spectroscopic quality redshifts. 

Fig. 1 shows the headline constraints on dark energy (left panel) and deviations from General 
Relativity (right panel). Table 1 gives details of the assumed area and number density of sources 
for the different surveys we forecast. 

The WGL surveys have photo-z quality redshifts. SKAl early is assumed to get these from 
a DES-like survey so we assume photo-z error, 5^ = 0.07(1 -|-z) and 5 tomographic bins of equal 
number density out to z = 2. SKAl and SKA2 are assumed to get redshifts from a Euclid-like 
survey with 5^ = 0.05(1 +z) and 10 tomographic bins of equal number density out to z = 2. 
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SKA1 5000 de^ with Planck 



SKA2 3pi with Planck 



SKA1 5000 de^ with Planck 



SKA2 3pi with Planck 



Figure 1: Constraints on dark energy [left panels] and deviations from GR [right panels] for SKAl over 
5,000deg^ [top panels] and SKA2 over 30,000deg^ [bottom panels] including Planck priors. Black el¬ 
lipses show photometric WGL constraints only. Red ellipses show spectroscopic LSS constraints only. 
Blue ellipses show the combination of WGL and LSS including cross-correlations. Cyan ellipses show this 
WGLxLSS constraint combined independently with an SKA intensity mapping (IM) survey using the MID 
instrument and the green contours the same but with the SUR instrument. All constraints are 68% confidence 
contours. 


The LSS surveys have galaxy redshift distributions described in the HI bias and simulations 
chapter of this book. Our forecasts assume 20 tomographic bins up to z = 0.6 for SKAl and 
40 tomographic bins up to z = 2.0 for SKA2. We use the exact C{1) formalism (not the limber 
approximation) and include the effects of Redshift Space Distortions (RSDs) according to the for¬ 
malism of Kaiser (1987). Both these effects are neglected in the WGL forecasts because the broad 
tomographic bins make their impact negligible. 

Our FM analysis forecasts constraints for a set of cosmological parameters: 
{flm,0^,0£)£,>vo,>Va,/j,a8,ni,(>,2o;2o(l +^o)/2}- As well as the standard wCDM parameters, 
is a free amplitude on galaxy bias for each shell and 2o) Go(l +7^o)/2 are parameterisations of 
deviations to General Relativity that modify the Poisson equation and the ratio of metric potentials, 
our ability to constrain these parameters quantifies our ability to test gravity on cosmic scales, see 
Kirk et al. (2013) for more detail. When quoting constraints on dark energy we marginalise over 
the cosmological parameters and galaxy bias but keep the modified gravify paramefers fixed. When 
quofing consfrainfs on modified gravify we marginalise over cosmology, including wq and Wa, and 
galaxy bias. Priors consisfenf wifh fhe lafesf Planck femperafure consfrainfs are included and fhe 
same used in fhe BAG chapfer in fhis science book. 

If is worfh remarking fhaf our forecasting approach is nof fhe only mefhod fhaf has been used 
fo estimate fhe join! power of a phofomefric WGL survey and a specfroscopic LSS survey. We 
freaf all our observables as projected angular power specfra, C(£)s. The alfernafive is fo model fhe 
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LSS survey with a full 3D P{k,z) analysis then combine with WGL C(^)s. Our projection over the 
line-of-sight width of a tomographic bin throws away information compared to a full 3D analysis. 
For a photo-z WGL survey the photo-z scatter and the broad lensing kernel combine to produce 
a redshift resolution floor, below which precision cannot be pushed. This means our relatively 
coarse redshift binning captures all the available information. The spectroscopic survey has no 
such limitation, this is why we use a very large number of narrow tomographic bins to capture 
equivalent redshift-evolution information as contained in a full 3D analysis. The advantage of our 
approach is that all observables are in the same form and cross-correlations arise naturally. Mixing 
a 2D and 3D analysis leads to some rather ad-hoc formalisms to combine different observables. 

Our results finding significant improvement from cross-correlation agree with Gaztanaga et 
al. (2013) who use a mixed 2D and 3D formalism. However, both Cai & Bernstein (2012) and 
de Putter, Dore & Takada (2013) adopt a similar mixed approach but see little improvement from 
cross-correlation. This is an active area of research and discussion on the correct formalism or even 
how to accurately compare results is continuing. 

DE Figure of Merit (FoM, see Albrecht et al. (2006)) for our various survey combinations are 
shown in Table 1. The joint WGFxFSS constraints are extremely strong. While the FSS surveys 
for both SKAl and SKA2 favour large areas over depth, the joint WGFxFSS constraint for SKA2 
prefers the deeper but smaller 5,000deg^ survey over the shallower 30,000 deg^ survey. The Modi¬ 
fied Gravity constraint shown in the right panel of Fig. 1 has a WGFxFSS constraint that is orders 
of magnitude stronger than either probe alone. This is due to pronounced degeneracy breaking from 
the combination of a probe using light (WGF, sensitive to the sum of metric potentials 'P-|-<I>) and 
a probe using galaxies as non-relativistic tracers (FSS, sensitive to the Newtonian potential 'P). 

SKAl WF -I- FSS represents a formidable dataset but it will of course not be the only one 
available. Dedicated optical surveys such as DES and eBOSS will be available and relatively ma¬ 
ture by the completion of SKAl. In terms of number density and area a representative SKAl 
survey configuration is comparable to DES-i-eBOSS with perhaps marginally better photo-z accu¬ 
racy achieved with SKAl. We find that SKAl (WGF-i-FSS-i-IM) provides comparable constraints 
on DE to a representative DES-teBOSS forecast while SKAl significantly outperforms (more than 
a factor of five improvement) DES-teBOSS in constraining deviations from GR, probably driven by 
wider sky coverage for the spec-z survey. The fact that both probes are collected by the same instru¬ 
ment offers benefits for the understanding and control of systematics. Of course having datasets of 
comparable statistical power from both optical and radio surveys is of considerable scientific value 
of itself. It is with SKA2 that the project becomes definitively world-leading with much higher WF 
source density than the Euclid mission and a spectroscopic quality FSS survey covering the full 
area and depth of the WF measurement. 

4. Cross correlations for redshift calibration. 

Photometric redshifts are less accurate but much faster and cheaper to gather than spectro¬ 
scopic redshifts. SKA will allow a large number of spectroscopic-quality redshifts to be gathered 
as it measures both galaxy position and redshift via the 21cm line for low redshift galaxies. For 
higher redshift galaxies used for example in the WGF survey it will require photometric-quality 
redshifts gathered by other surveys like Euclid. 
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A photometric redshift from a Euclid-like experiment can have a random scatter compared 
to the true redshift of order 8^ ~ 0.05(l-i-z). As well as this Gaussian scatter on the true redshift 
(which can lead to objects being assigned to the incorrect tomographic bin). We note that the true 
redshift distribution can become much broader than the extent of the corresponding tomographic 
bin. There are problems estimating the true mean redshift of a certain tomographic bin and what 
are known as “catastrophic outliers”, galaxies whose redshift has been severely misidentified due 
to failures in the photo-z estimation pipeline, usually these are defined as estimates more than 3a 
away from the true redshift, which e.g. can occur at optical wavelengths if the the Lyman and 
Balmer breaks are confused. Systematic shifts in the redshift distribution of tomographic bins can 
induce significant biases in cosmological analysis, so that the mean of these distributions needs to 
be known to better than a few parts in a thousand. 

Cross-correlation with a spectroscopic survey which covers some or all of the same galaxy 
distribution as the photometric survey can be used to identify and mitigate these errors and calibrate 
the photometric redshift distribution. 

We repeat the WGLxLSS forecasts made above but include now a range of nuisance param¬ 
eters, {5,.,, Af?^which quantify our uncertainty on the photometric redshifts. Here & , is the 
Gaussian uncertainty on the photo-z estimate in each photo-z bin, i, with fiducial values 0.05 and 

is the bias on the mean redshift of bin i due to photo-z mis-estimation with mean values 0. We 
ignore catastrophic outliers in this analysis. Each of these parameters is allowed to vary in our EM 
analysis and the inclusion of the WGLxLSS cross-correlation allows them to be constrained. 

See Eig 2 for constraints on these nuisance parameters for the SKA2 371 survey. Cross¬ 
correlation of the WGL survey with the spectroscopic ESS survey improves our estimate of the 
photo-z distribution by up to factor of 10. We emphasise that this is a cross-correlation between 
photometric WGL and spectroscopic ESS, in future we intend to study the cross-correlation be¬ 
tween a ESS analysis using the photometric WGL galaxies and the spectroscopic ESS survey. This 
should prove even more effective at calibrating the photo-z errors because it is not hindered by the 
broad geometric kernels which restrict redshift resolution in WGL. 

5. Cross-correlations with CMB lensing 

The gravitational potential of the large scale structure generates a deflection of the trajecto¬ 
ries of the Cosmic Microwave Background (CMB) photons. This effect, known as CMB lensing 
consists of a remapping of the underlying unlensed temperature and polarization fields. Mafhemat- 
ically, CMB lensing is described as follows: We introduce a vector field d(n) (the deflection field) 
such that the lensed temperature T (n) and unlensed temperature f (n) are related by 


r(n) = f(n + d(n)) 


(5.1) 


and analogously for the Stokes parameters 2(n), U (n) which describe linear CMB polarization. 

To lowest order in perturbation theory, the deflection field d(n) is the gradient of a scalar 
lensing potential (i.e. d(n) = V0(n)) which can be written as a line-of-sight integral: 



(5.2) 
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SKA2 3pi 



SKA2 3pi 



Figure 2: Constraints on photometric redshift error [top] and bias on mean redshift [lower] as a function 
of mean redshift of each tomographic bin for SKA2 over 30,000deg^. WGL- only constraints are shown in 
black, WGLxLSS constraints in blue. Constraints including Planck priors are shown as dashed lines. 


where is the Newtonian potential, rj is conformal time, T]rec is the epoch of last scattering, and 
X is the angular diameter distance in comoving coordinates, see Lewis & Challinor (2006) for a 
review). 

CMB lensing modifies the Gaussian structure of the primary anisotropies and generates a 
correlation between the temperature and its gradient Hu (2000). These couplings can efficiently be 
used to construct an estimator, quadratic in the observed temperature, that can be applied to data to 
recover the lensing potential Okamoto & Hu (2003). 

High resolution and high sensitivity CMB experiments can therefore provide a new cosmolog¬ 
ical probe of the large-scale structure of the Universe that is complementary to that obtained from 
galaxy surveys. Indeed CMB lensing is mostly sensitive to structure located in the 1 < z < 5 range. 

The lensing potential reconstructed from CMB lensing can be thought of as the projection on 
the sky of all the mass distribution up to the last scattering surface. As such we expect significant 
angular cross-correlation with the large-scale structure observables of radio surveys. The redshift 
dependence of the CMB lensing kernel means that, to ensure a strong cross-correlation, we ideally 
want a LSS survey with many galaxies at z 3> 1 and a WGL survey with many galaxies at Z 2. 
SKA will produce very useful surveys in this respect, more useful than Euclid for example when 
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Survey 

Area [deg^] 
Photo-z 

Planck 

WGL LSS 

SPT3G 

WGL LSS 

SKAl early 

5,000 

5.7% 

4.1% 

2,8% 

1.7% 

SKAl 

5,000 

4.8% 

3.4% 

2.2% 

1.3% 

SKAl 

371 

3.8% 

2.2% 

4.3% 

1.9% 

SKA2 

5,000 

2.9% 

2.4% 

0.9% 

0.7% 

SKA2 

371 

1.6% 

1.3% 

1.1% 

0.9% 


Table 2: Constraints on the amplitude of the cross-correlation between the Planck and SPT3G lensing 
potential and the SKA surveys observables: weak lensing (WGL) and galaxy clustering (LSS). 

considering this cross-correlation with CMB lensing. 

The first reported detection of the gravitational lensing of the Cosmic Microwave Background 
was made by correlating WMAP data with the radio galaxy counts from the NRAO VLA sky survey 
(NVSS) (Smith, Zahn & Dore (2007), Hirata et al. (2007)). With the advent of arcminute scale 
CMB experiments (ACT, SPT), and full-sky CMB data from Planck, prospects for cross-correlation 
between CMB lensing and future radio data are extremely encouraging. 

SKA will provide both precise measurement of the position and the shape of galaxies. In the 
following we use the same redshift distributions assumed in the Weak lensing chapter from this 
science book, and consider a value constant for the bias b = \. 

On the CMB lensing side, we consider the following surveys: 

• Planck. We consider the Planck lensing potential from the Planck 2013 results (Planck col¬ 
laboration (2013)). The current Planck lensing map should be replaced in October 2014 and 
will include the full Planck data. We thus multiply the 2013 lensing noise levels by a factor 
0.8 to account for this added data. 

• South Pole Telescope (SPT). The SPT collaboration will observe the CMB polarization 
anisotropies to arcminute resolution on a 2500 sq. deg. patch in the southern hemisphere. 
The resulting lensing map will then be of a much better quality than the Planck one, but on 
a smaller area on the sky. The forecasted noise levels were provided by G.Holder and G. 
Simard (private communication). 

In the following we quantify the level of detection of the cross-correlation between CMB lens¬ 
ing and the SKA observables. We use the values indicated in table 1 for the different incarnations 
of the SKA surveys. Since we are considering cross-correlations with external data set, we need to 
restrict the sky fraction to the maximum common area between SKA surveys and Planck/SPT data. 
We therefore restrict the SKA surveys to 2500 sq. deq. when correlating with the SPT lensing map. 
For the correlation between SKA surveys and the Planck lensing map we consider a maximum area 
of 15000 sq. deg. 

Constraints on the amplitude of the cross-correlation between the SKA surveys observables 
(weak-lensing and galaxy clustering) and the CMB lensing potential are summarized in table 2. 
Some examples of those cross-correlations are shown in fig. 3. 


9 







Cross correlation surveys with the SKA 


Filipe B. Abdalla 



500 1000 1500 2000 


e 

Shear - CMB lensing 



■■ 5KA2_3pi • Planck 
SKA2_3pi - 5PT3G 



3.0 



500 1000 1500 2000 


e 



500 1000 1500 


£ 



Figure 3: Weak-lensing (left column) and galaxy clustering (right column) for SKAl-5000 (top row) and 
SKA2-3;r (bottom row) cross-correlations with the CMB lensing potential. Errors bars correspond to the 
Planck (green) and SPT3G (blue) noise levels. 


Prospects for cross-correlations between SKA and CMB lensing are very encouraging. As 
can be seen from the plots in fig. 3, most of the signal will come from the use of high resolution, 
small scale CMB data such as SPT3G (similar results would hold using the Atacama Cosmology 
Telescope (ACT) specifications). This can be seen as the blue errors bars being much smaller 
than the green ones, which represent the correlation with Planck. However, thanks to its large 
sky coverage, Planck will dominate on the largest scale, despite being intrinsically noisier than 
arcminute scale CMB data. 

The cross-correlation of SKA probes with CMB-lensing will improve our control of system- 
atics and our ability to measure cosmology. CMB-lensing is sensitive to the integrated matter 
density between the observer and last-scattering without the mediation of biased tracers, helping 
us to control galaxy bias in a similar way to WGL but because the measurement of CMB-lensing 
is independent it can also calibrate important WGL systematics including shear measurement bias 
and intrinsic alignments, see Vallinotto (2013) for an example of the kinds of improvements that 
are possible. 
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6. Conclusions. 

In this chapter we have investigated the synergies from cross-correlating different SKA datasets 
as well as cross-correlation with CMB lensing datasets from other missions. We conclude that, us¬ 
ing internal cross-correlation, the SKA will be able to calibrate the redshift inaccuracies present in 
its weak lensing sample, optimising the possible statistical measurements of cosmology. We see 
from our results that LSS measurements with galaxies from the SKA are not very competitive dur¬ 
ing phase 1 but they are still important as they help with the aforementioned redshift calibration, 
particularly the planned ~100 deg^ deep survey which will overlap fully in redshift with the WGL 
survey source population. We estimate this calibration to be possible at the sub-per cent level with 
SKA2. 

We have shown that these cross-correlations provide huge gains for our Dark Energy anal¬ 
ysis with SKA phase 2, but the largest gains come when we study Modified Gravity where, for 
SKAl and SKA2, gains of several orders of magnitude are possible from combining WGL and 
LSS datasets. 

We also present constraints on the amplitude of the cross correlations between the lensing sig¬ 
nal of future CMB experiments and the SKA and a sub present level constraint on the amplitude 
is possible. This would impose stringent constraints on the bias of SKA galaxy samples, calibrate 
systematic uncertainties in our WGL measurement and help constrain other cosmological parame¬ 
ters. 
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